function [ blip, tout ] = get_blip_from_fq(F, amiga_model, ... 
    filter_method, length)

filter_short_len = get_filter_short_len();
fs = get_sampling_fq();

if strcmpi(filter_method, 'fq')
    F = filter_fq(F);
end

blip = real(ifft(F));

if strcmpi(length, 'short')
    tout = linspace(0, (filter_short_len-1) / fs, filter_short_len);
    blip = blip(1:1:filter_short_len);
elseif strcmpi(length, 'full')
    tout = linspace(0, (fs-2)/fs, fs-1);
end

if ~strcmpi(filter_method, 'fq')
    blip = fir_min_phase(blip);
end

if strcmpi(filter_method, 'conv')
    blip = filter_impulse(amiga_model, blip, filter_method);
end

blip = blip / trapz(blip);

end